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Abstract 

We present an exact nonperturbative flow equation for the average action for quarks 
which incorporates the effects of gluon fluctuations. With suitable truncations this 
allows one to compute effective multiquark interactions in dependence on an infrared 
scale k. Our method amounts to integrating out the gluons with momenta larger 
than k. 



1 Introduction 



Effective multiquark interactions have often served as a starting point for QCD- 
motivated investigations of hadronic physics. A well-known example is the Nambu- 
Jona-Lasinio model which is based on an effective four-quark interaction. The 
idea is that these models are valid for quark fluctuations with momenta smaller than 
a given scale k whereas fluctuations with higher momenta and gluons are thought 
to be already integrated out. In the more fundamental context of QCD one would 
like to have some information on the strength of such an interaction in dependence 
on k. Furthermore, in a realistic theory the effective four-quark interactions will 
be momentum-dependent. For heavy quarks this momentum dependence translates 
into the form of the heavy quark potential which is the basis for the non-relativistic 
quark model for charm and bottom-bound states 0. Recently a QCD- motivated 
effective four-quark interaction (at a scale k ^ 1.5 GeV) has also been used to 
derive the emergence of mesonic bound states and chiral symmetry breaking from 
the solution of nonperturbative flow equations 0]. In order to make direct contact 
with QCD all these approaches need a method capable of integrating out the gluonic 
degrees of freedom. 

The problem of integrating out completely the gluons seems almost as difficult as 
the full solution of QCD. Even if an effective infrared cutoff is present in the quark 
sector the gluons will lead to strong nonlocalities in the quark interactions. Formally, 
one would have to derive the effective action for gluons and quarks, if needed with 
an additional infrared cutoff ~ /c in the quark sector. From there the gluonic field 
equations can be derived, giving the (classical) gluon field A^^''^ as a functional of 
the quark field ip. Inserting this solution into the effective action then leads to an 
effective action for the quarks from which the effective multiquark interactions can 
be extracted. Obviously such a program must become very complicated since due 
to confinement a purely gluonic description of the gauge degrees of freedom seems 
not very appropriate in the low momentum range. New non-perturbative effective 
degrees of freedom are expected to become relevant for very low momenta below the 
confinement scale A. 

For many practical purposes, however, the extreme low momentum behaviour 
of the effective quark interactions is actually not needed. For example, the use of 
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the effective heavy quark potential for the low lying charmonium and bottonium 
states essentially needs information about the effective four-quark interaction at 
momenta ~ (300 MeV)^. Since the quark momenta act as an effective infrared 
cutoff also for the gluon fluctuation^, we do not expect that confinement effects 
are dominant in the range of sufficiently large momenta. One may therefore think 
of using there an approximation with an additional infrared cutoff k both for the 
gluon and the quark degrees of freedom. This idea is realized by the concept of the 
average action which obtains from integrating out all quantum fluctuations 

with momenta larger than k. For a given A; > one does not expect that always 
describes very accurately the n-point functions for momenta smaller than k - m 
particular the confinement physics in QCD is cut out for large enough k. On the 
other hand the effective average action can be interpreted as a good approximation 
to the generating functional for IPI Green functions with momenta sufficiently larger 
than k. A possible first approach for a computation of the potential in the range 
relevant for heavy quarkonia could therefore be attempted by a computation of the 
quark average action F^. [■?/'] for k somewhat larger than the confinement scale. We 
will see next how this can be improved by lowering /c — 0, provided the momenta 
in the effective quark n-point functions are sufficiently large. 

The dependence of the average action on the scale k is described by an exact 
non-perturbative flow equation |p. This functional differential equation needs a 
truncation in order to become approximately solvable in a non-perturbative context. 
Typically, one may retain only gluonic two-, three- and four-point functions in the 
gauge sector. The question arises in what context this type of truncation can lead to 
meaningful results. First we should point out that such a truncation is not restricted 
to the range of applicability of perturbation theory. It is general enough to account 
for a nonperturbative form of the gluon propagator and vertices. Nevertheless, 
a simple gluonic description is not expected to remain valid once effective non- 
perturbative degrees of freedom different from the gluons play an important role. 
(This is partly related to situations where higher n-point functions for the gluons 
become important.) The range in k where a "gluonic truncation" of the average 
action can be trusted is limited by this fact. For k larger than the confinement scale 
the precise description of the confinement mechanism should not be of dominant 

^We work here in the Euchdean regime where > and quarks are ofF-shell. 
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importance. For k as large as a few GeV one expects that a perturbative description 
of Ffc becomes valid. It is well conceivable that once k is lowered, there first appears a 
region where perturbation theory breaks down but a truncation to a few gluon Green 
functions still remains a valid approximation. Let us now discuss what happens if one 
tries to follow the flow of the average action beyond this region towards — > within 
such a "gluonic truncation" . One expects that the low momentum behaviour of the 
n-point functions becomes only poorly described and may even give a qualitatively 
wrong picture. In contrast, the Green functions with momenta ~ (300 MeV)^ 
should only mildly be affected by the error from the insufficient truncation. Since 
an infrared cutoff is already present due to the nonvanishing momenta the variation 
of an additional and smaller infrared cutoff k should not be very relevant anymore. 
We conclude that a "gluon truncation" can give reliable results even for A; — > (all 
fluctuations included), provided the momenta of the Green functions are sufficiently 
large. In particular, one wonders if a computation of the heavy quark potential 
V{R) for the range in R relevant for the low-lying bound states may be feasible 
along these lines. 

Let us next turn to the light mesons. The scale for the formation of mesonic 
bound states was found to be around k^ ^ 700 MeV. For k < k^ a. mesonic 
description becomes appropriate [Q]. Following the flow equation of an effective 
quark-meson model for k < k^ one finds that chiral symmetry breaking sets in 
at A; ~ 450 MeV. Subsequently the quarks acquire a constituent mass ^ 350 
MeV from their Yukawa coupling to the chiral symmetry-breaking order parameter. 
This implies that quarks effectively decouple for k ~ 300 MeV. Since mesons do 
not carry colour one may hope that the details of the confinement mechanism are 
not crucial for an understanding of the meson mass spectrum and decays. This 
may hold despite the fact that already the formation of light meson bound states is 
certainly a highly non-perturbative effect. The exchange of many gluons needs to 
be considered for the description of bound states. 

In order to understand this issue more precisely one may imagine for a moment 
that at some high scale ki the classical gluon exchange gives rise to an "initial" 

— 2 — 

effective four-quark interaction ~ tp'1^4'^'4''ltJ.''P ^ind the gluons are neglected for 
k < ki. If one follows the flow equations for the fermionic model for k < k^ the 
presence of an "initial" four-fermion coupling influences the trajectories and results 
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for low enough A; in a complicated structure of the effective four-quark interaction. 
In this sense a great part of "multi-gluon exchanges" is already incorporated in the 
solution of the purely fermionic flow equation. Nonperturbative gluon effects arise 
even though the gluons do not appear explicitly for k < ki. This is, however, not 
the whole story: There are also contributions from gluon fluctuations which are not 
included in the standard flow of an effective fermionic theory. Contributions from 
those "residual gluon fluctuations" lead to important corrections to the standard 
fermionic flow equations. They will be computed explicitly in the present paper. 
Corrections to the flow of the four-quark interaction involve the two-, three- and four- 
point functions for gluons. In this context the truncation in the gluon sector only 
enters through the "residual gluon contributions". In an appropriate formulation 
the "residual gluons" influence only very indirectly the mesons via their effects on 
the quarks. For k < 300 MeV their contributions are suppressed by the quark 
constituent masses. It is in this type of formulation that the precise description of 
conflnement may not be needed, and even a relatively inaccurate "gluon truncation" 
may not destroy the reliability of a calculation of meson properties for A; ^ 0. 

In order to address the questions raised in this introduction more quantitatively 
one needs a method which integrates out the gluon fluctuations with momenta 

> k'^. Formally, it is easy to do so at a given scale ki. One needs to compute 
the effective average action for quarks and gluons F^J^jA] at this scale. Solving 
the classical field equation for the gluon field A in dependence on 'ip and inserting 
this classical solution into rfci[-(/^, A] yields exactly an effective average action FfcJ-?/'] 
involving only the quark fields. The quark-effective action F^^ [ip] may then be used 
as an initial value for solving a pure fermionic flow equation for k < ki. Obviously, 
the shortcoming of such an approach ist the complete omission of the effects of 
residual gluon fluctuations with < kf. Choosing a different scale k2 for eliminat- 
ing the gluons will lead to a different result. Such a sharp transition between the 
quark-gluon system and a description involving only quarks necessarily introduces 
a certain degree of arbritariness. 

We propose here a more reflned method which changes the classical solution 
for the gluon fleld in the course of the evolution towards lower k, thus reflecting 
the change in the form of FfefV',^!]- The result is a smooth procedure for integrat- 
ing out the gluons, where at every scale k all gluon fluctuations with > k'^ are 
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included. Nevertheless, we obtain a flow equation for the quark effective average 
action [■?/'], where gluon flelds do not appear explicitly. The inclusion of contribu- 
tions from residual gluon fluctuations induces correction terms in the flow equation 
for the quark interactions. In particular, we choose here a formulation where the 
/c-dependent classical solution for A as a functional of includes an effective infrared 
cutoff ~ /c. As a consequence, the only nonlocalities in Tk[ip] concern length scales 
shorter than k~^. For the fermionic low momentum modes Ffc[^] is an effectively 
local action. For example, a derivative expansion is meaningful for <^ k^. The 
expected nonlocalities arising from the complete elimination of gluon fields (for ex- 
ample a four-quark interaction ~ ^ appearing already in the Born approximation) 
build up only step by step as A; is lowered to zero. As a result of this method we 
will end with an exact nonperturbative fiow equation for the scale dependence of 
Ffef'^]. Approximations will be needed to solve this equation but they are not lim- 
ited to perturbative concepts. For example, we may neglect interactions which 
should be roughly equivalent to neglecting the influence of the baryons on the me- 
son properties. Also the correction terms reflecting the residual gluon fluctuations 
require limited knowledge about the effective gluon propagator and vertices. It is 
hoped that rather crude approximations for the gluonic vertices can already lead to 
satisfactory results. 

In sect. 2 we develop our formalism for a simple model of two scalar flelds. The 

flow of the effective average action for one of the scalar flelds obtains by integrating 
out the other scalar fleld at any scale k. This is generalized to quarks and gluons in 
sect. 3. In sect. 4 we give a flrst demonstration how this formalism describes the 
flow of the two- and four-point function in the effective quark theory. The special 
case of heavy quarks is addressed in sect. 5. Here we argue that only the evolution of 
the gauge fleld propagator is needed in this limit and we compute the corresponding 
flow equation in one of the appendices. Finally, our conclusions are contained in 
sect. 6. 

2 Reduction of degrees of freedom 

In this section we consider for simplicity two types of scalar flelds, (p and ip. We want 
to develop a formalism how to translate evolution equations for the effective average 
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action for Lp and i}) into corresponding equations involving only The reader may 
associate ip with the gluon fields and ip with the quark fields. Our aim is then the 
construction of the effective average action for quarks out of the coupled quark-gluon 
system. This amounts to integrating out the gluonic degrees of freedom represented 
in the simplified model of this section by ip. We start with the scale-dependent 
generating functional for the connected Green functions 

Wk[J, K] = \nJ D^'D^' exp -{S[<^', + /^^^ S[v'] + A?^5[^'] - JV - ^V} 

(2.1) 

Here we denote the degrees of freedom contained in ip' (for example the Fourier 
modes) by y^'" and similar for ip' , J and K, withQ 

jt<^' = kH' = k;^'^ (2.2) 

We have introduced an infrared cutoff quadratic in the fields 

^'i^^SW] = (2.3) 



and similar for tp. This suppresses the contribution of fiuctuations with small mo- 
menta < k"^ to the functional integral (|2.1| ). Typically R^^\r^^^ are functions of 
with the properties 

, fc^ for <^ /c^ , ^ 

^^-^ 2 , ,^ . 2 ,2 (2.4) 

g exp(-|2) for k'^ 



as, for example, 



Rl^> = — 1^4^ (2.5) 

l-exp(-g) 

The effective average action Tk[ip, ip] is related to the Legendre transform of K] 

tk[ip,tp] = -Wk[J,K] + J^ip + K^ij (2.6) 

^" = ^ iP^=^-^ (27) 



by subtracting the infrared cutoff term 



r,[^,,p] = r,[^,,p] - A'f^Slip] - A^SliP] (2.8) 



^We use indices a, a' etc. for ip and /3, (3' etc. for ip. 
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For k the infrared cutoff A^S = A^j^^ S + A^ff^ S vanishes and Tq is the usual 
generating function for the IPI Green functions. Using the quadratic form of AkS 
it is straightforward to derive an exact non-perturbative evolution equation for the 
dependence of the effective average action on the scale k {t = Ink) 



:2.9) 



Here f^^^ = if ^ + Rk and the inverse propagator r^.^"* is the second functional 
derivative of with respect to the fields. The matrix Rk = R^^^ + -R[!^^ is block 

~ (2) 

diagonal in (p and i/j spaces. The presence of the infrared cutoff Rk in guarantees 
infrared finiteness for the momentum integral implied by the trace even in case of 
massless modes. Ultraviolet finiteness is guaranteed by the exponential decay of 
dRk/dt ( |2.5| ). A solution of the flow equation (|2.9| ) interpolates between the classical 
action for k oo (or k equal to some ultraviolet cutoff A) and the effective action 
for k —>■ 0. 

The generating functional for the connected Green functions for obtains from 



(13) for J = 

Wk[K]=Wk[J = 0,K] (2.10) 
Correspondingly, we may introduce an effective action expressed only in terms of ip 

fkm = fk[ipkm,^] (2.11) 

= Tk[Mi^],^] + A^^^S[M^]] (2.12) 
by inserting the fc-dependent solution of the field equation 

= (2.13) 

This defines ipk as a fc-dependent functional of ip. Since for every scale k the definition 
( p. 61) implies 

and therefore, for J = 0, 



dipP dipP di/ji^ dip'' 1^=^^ 



+ ZT Z: ' = K; (2.15) 
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it is easy to verify that Tk[ip] is the Legendre transform of ^^^[-ft'] ( |2.1(j| ). One 
concludes for — that Tq[iIj] is the generating functional for the IPI Green 
functions for tp. 

We want to employ the flow equation ( p.9| ) for finding the /c-dependence of [■?/;]. 
In addition to the corresponding equation for only one type of fields we have here 
additional contributions from the fc-dependence of A^j^^ S in (|2.1|) . We prefer to keep0 
this additional infrared cutoff for the "gluons" in order to maintain approximate 
locality of T^lip] on length scales large compared to k~^. The evolution equation for 
Ffc [ip] can now be obtained from a variable transformation which amounts to a shift 
of (f around ipk[ip] 

r = ^''-^tm (2.16) 

Using the general formalism of appendix A one finds 




It is easy to verify that this equation reduces in the limit = to the equivalent 
of eq. ( ^.91 ) for fields ip only. The corrections in the first two terms involve the 
explicit form of the "classical solution" v^a:[V']. If one is interested in IPI Green 
functions for ip with a given number of external legs one only needs a polynomial 
expansion of (fk [4^] up to a given order. For example, the evolution of the term ~ ip'^ 
in Tk[ip] needs the classical solution up to the order ip^ if the series ipk[i'] starts 
with a term quadratic in ip. Additional knowledge of the form of Ffc[(^,^/'] beyond 
its value for = is only needed for the last correction term in the form of 

Only the ^'-dependence of the effective (^-propagator plays a role for the study of 

IPI functions for ip. 

different approach uses separate cutoff scales k and k for the "gluons" ip and "quarks" 
ip. Letting k ^ for fixed k one completely integrates out the gluons. This alternative method is 
most appropriate for massive fields. Then the mass term acts as an infrared cutoff and prevents 
the appearance of strong nonlocalities for < k"^ 
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We finally note that the definition of [■?/'] proposed in this paper is not the only 
possible choice of an effective fermionic action. An alternative formulation is briefly 
discussed in appendix B. 



3 Integrating out the gauge fields 

The generalizaton of the discussion of the last section for quarks and gluons follows 
ref. 1^. All formulae of the last section remain valid for arbitrary bosonic flelds 
if the indices a, [3 ( |2.2| ) include internal indices and Lorentz indices in addition to 
momentum labels^. If is a Grassmann variable as appropriate for fermions the 
matrix Rk in (|2.9|) becomes (cf. appendix C) 

tik — rik - I'^.ij 

Also ip* should be replaced by if) and the index summation over (3 should involve 
both ip and ip separately.^ For the gauge flelds we will choose here a formulation with 
explicit ghost variables in close analogy, but slightly different from the formulation in 
ref. 0. This makes our formulation as close as possible to the language of standard 
perturbation theory. 

We start with the action including a gauge-flxing term in the background gauge 
and a corresponding action for the anticommuting ghost flelds ^, ^ 

SW, i\ a; A] = S[iP', A'] + S,f[a; A] + S,hW, a; A] (3.2) 

Here S* is a gauge invariant functional of the fermion flelds ip, ip and the gauge fleld 

A'^ = A^ + a,. (3.3) 

The background gauge fleld A^ appears in the gauge flxing and ghost terms 

1 r,. 



Sgf = —jd'xGlG^ (3.4) 
= m^Wyal (3.5) 
Sgh = Jd'xey{-D>^[A]D^[A + a])y,e' (3.6) 



'^For a complex field (p the sums are also over negative internal indices i with ip '((?) = (fl {q) ■ 
^One may again use negative internal indices for labeling ip, i.e. ip~^{q) ~ ^Ji{q) and use the 
summations of the preceding section including the negative indices. 



9 



Here -D^[^] is the covariant derivative in the adjoint representation in presence of 
the background gauge field A. The generating functional for the connected Green 
functions is defined as usual 

W[7], C, K;A]= J Vij'ViP'V^'V^'Vaexp -{S - J d'^x[f]%l)' + r/V^' + C^' + + Ka]} 

(3.7) 

We have introduced here also sources C for the ghost fields and note that the source 
couples to the gauge field fiuctuation and therefore transforms homogeneously 
under gauge transformations as an adjoint tensor. The /c-dependent version Wk 
obtains from W by adding to S the infrared cutoff piece 

^kS = A?)^ + Ai^)^ + aI^'^^S (3.8) 

Here the fermionic cutoff reads 



1) / d'x^l,'Z^^Sl'D,[A]) r''^\-D'[A]/eW (3.9) 

with the covariant derivative in the appropriate representation {D'^ = D^D'^) and 
a dimensionless function. The factor 1/2 is appropriate for Majorana fermions 
0,0. For the gauge field cutoff we choose 



V[A] r 



{A) 



al (3.10) 



vz 



with T>\A\ an appropriate operator generalizing a covariant Laplacian in the adjoint 
representation which will be explained below. The matrix ^a./c accounts for an 
appropriate wave function renormalization. Finally, we take for the ghosts 

Af^-s = e;,(i?f))^;r 

= / d'^x%\Z,^^kVAA\Tt\v,\A\lk^)\\e (3.11) 

with Ps[v4] = — Z}^[y4] in the adjoint representation. A good choice for the dimen- 
sionless function isQ 

rAv) = (3.12) 

1 — e y 

^The function r'^'^ may also be chosen differently from (3.12) in order to avoid that Rk diverges 
for vanishing covariant momenta. 
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such that 

hm Rk = Zkk'^ (3.13) 

The fc-dependent functions Z^^k,2A,k and Zgh^k will be adapted to corresponding 
wave function renormalization constants in the kinetic terms for the fermions, gauge 
fields and ghosts. In principle, they can depend on the background field A. The 
infrared cutoff piece AkS cuts off all quantum fluctuations with covariant momenta 
smaller than k in the functional integral defining Wk- For covariant momenta larger 
than k the infrared cutoff is ineffective and its contribution to the propagator is 
exponentially suppressed. 

Performing a Legendre transform and subtracting the IR cutoff piece again (c.f. 
(p.6[), (|2.8|)) we arrive at the effective average action Tklip, ^, A, A], where A = A + a 
and a is conjugate to K. The dependence of on the scale k is described by 
an exact evolution equation analogous to eq. ( p.9|) , with a negative sign for the 



contributions ~ R^^^ and -R^f^''. We note that F^ only involves terms with an even 
number of ghost fields due to the symmetry ^' — > — f', ^' — > of the Sgh and 
A^^'^^S'. In consequence, the ghost field equations 

^ = 0.^ = (3.14) 

have always the solution ^ = C, =0. We therefore can extract the propagators and 
vertices for the physical particles from the effective action for ^ = ^ = 0: 

Tk[ip,A,A] = Tk[tp,0,A,A] (3.15) 

Nevertheless, the evolution equation for Tf^lip, A, A] obtains a contribution from the 
variation of the infrared cutoff of the ghost fields as given by 

-(i)Tv{(|fl?')(rf + fi.)")-<^. (3.16) 

iv|(|i?r)(rr<'' + 4"'")"} (3.17) 



Here v'-^^+Rk in ( |3.16| ) is the matrix of second functional derivatives of Ffc + A^f ■'5' + 
A^.'^-'^' with respect to and A at fixed A. For this we have exploited that the larger 
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matrix of second functional derivatives of (Ffc + AkS)[ip,^,A,A] is block diagonal 
in the (ipyA) and ^ components for ^ = 0. The ghost dependence of rfc[-?/', ^, A, A] 
appears in the evolution equation for [■?/', A, A] only through the term Ek which 
involves the second functional derivative with respect to the ghost fields r^^^'"-^^ 
which is evaluated at ,^ = ^ = and may depend on ip, A, A.Q We will not pay much 
attention to the detailed form of r^,^^'"'^'* in the present paper and approximate it by 



its "classical" value (cf. ( p.6| )) 

^i9h)i2) ^ _jj^^A]D^[A]. (3.18) 

In order to complete the formal setup of our investigation we need to specify the 
operator V in eq. ( p.lO|) . A good choice is 

V[A] = ri'^^^^^A] (3.19) 

where F^^^*'^'' is the second functional derivative of Tk[ip, A, A] with respect to A for 
fixed A and ip = 0, evaluated at the point A = A. As in previous formulations [0, the 
effective average action [■?/', A, A] is gauge invariant with respect to simultaneous 
gauge transformations of ip, A and A. 

We can now apply the formalism of the last section in order to "integrate out" the 
gluon fields A. Derivatives with respect to if in the preceding section are replaced by 
derivatives with respect to A at fixed background field A. In particular, the classical 
field equation, whose solution is A^, reads (cf. (|2.13|) ) 



,1' ^ = 3.20 

At this point Ak becomes a functional of and A. Due to the covariance of the 
field equation ( 3.20 ) we may gauge transform^ any given solution into a 



corresponding solution [■?/'', A'] = A'^. In consequence, the effective action 

rj^jj,, A] = F4^, A = AlA]+ Aj^^SjA = Au[iJ,AlA] (3.21) 

^In the formulation of rcf. [Q the quantity r^^'*'"^^'' is replaced by For the computation 

of the evolution of Tk[4^, A, A\ performed in ref. @] this is equivalent to ( 3.18 ). 



^The functional dependence of Ak on ^ and A is such that a gauge transformation on -0 and A 
results in a corresponding inhomogeneous gauge transformation of A^. 
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is gauge- invariant. For our purposes we want to work with an effective action involv- 
ing only the quark fields. This requires to fix the background field A conveniently. 
A possible choice is 

A = 0, Ak[^P]=Af,[tlj,A = 0] (3.22) 
Tkl^] = r/^[^, A = 0] = Tk[tlj,A = Ak[^],A = 0] + Ai^'^SiA = A^^], A = 0]. 

In this version R^if^ becomes a simple function of momenta. Summarizing our adap- 



tation of eq. ( |2.17| ) for quarks and gluons and using the results of appendix C one 
obtains^ 



at- = - \-2j - "'I ^ ' ' ^'\ \ m ] ^ m>,. \ dt ] 



{AY 



+^(rfW..>i = >i*,U = oi + flr')-",(^ 

+ 2^^"'Hrl At-e,[i;,A,[i;],A = 0] (3.23) 



The factor (|) in the first term is absent for Dirac spinors. The flow equation ( p.23|) 
is the central equation of this paper. 

('2) 

In order to exploit this equation we need Aa;['?/'], [■?/', A^, 0] and [■?/', A^., 0]. 
Following ref. it is illustrating to divide Tklip, A, A] into a gauge invariant effective 
action depending only on one gauge field 

Tk[iJ,A]=Tk[ip,A,A = A\ (3.24) 

and a generalized gauge fixing term which vanishes for A = A 

Fr'^[V^,A^] = Tk[iJ,A,A]-Tk[^,A]. (3.25) 

Besides the classical gauge-fixing term and quantum corrections to it the part rf'^"^'' 
contains also the fc-dependent counterterms, as for example a fc-dependent gluon 
mass. (In perturbation theory this vanishes for k 0.) These counterterms are 



^For general A there is also a purely A dependent contribution to rfe['0, A, A] which arises, for 
example, from the effective infrared cutoff for the ghost fields [Q. This plays no role for the choice 
A = 0. 
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related to the gauge-invariant kernel [■?/', A] by generalized Slavnov- Taylor identi- 



ties 0, [T^. We will not discuss the role of the counterterms very explicitly in the 
present paper. (It is conceivable that at some later stage it may become advanta- 
geous to subtract them for an improved definition of [■?/'].) Using the explicit form 
of the gauge invariant infrared cutoff |0 

Ai^^^ = IJ dMA; - Al){R^)Z{Al - A^) (3.26) 

the relevant classical field equation (|3.2CI|) for A becomes 

+ (fi<-.)».(^j - + = 0. (3^27) 

As a consequence of the gauge invariance of Ffc['?/^,y4] the gauge degrees of freedom 
contained in A only get fixed by the last two terms in eq. ( p.27| ). Similarly, the 



inverse propagator for the gauge degrees of freedom obtains contributions only from 

^ k 

f1^) A, A] = f1^) [^,A]+ Fr^^(^) A, A] (3.28) 

(2) 

since the contribution from F^ '[iIj, A] vanishes for them. 



4 Effective quark interactions 



The flow equation (|3.23| ) is an exact nonperturbative equation. Its solution for k ^ 



contains the full information on all IPI Green functions for the fields tp. Its solution, 
however, will only be approximative and involves a truncation of the general form of 
Ffc['?/', A, A]. For a first illustration of the effective action for quarks F^, [■?/'] we make 
here the simple ansatz (with ip representing Dirac spinors) 



Ffc[^, A,A = 0] = I d^x {iZ^^ui'^^b^i, 

+ ^ZF,kFrF^^ + ^-Z^Ad.^n' + am} (4.1) 



where Ck[ip] is independent of A. Here F^^ is the field strength of a nonabelian 
SU{Nc) gauge theory and quarks are in the fundamental representation, with 

D.iP = d,,p-igkZ]!lAlT,^ (4.2) 
F;, = d,Al-d.Al + guZ]llf^/A;Al (4.3) 
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The k dependence of is encoded in the fc-dependence of the renormahzed non- 
abehan gauge couphng gk, the wave function renormahzations Z^^k, Zp^k and Ck[ip]. 
We will discard the fc-dependence of the ghost-action and take Ti^gh)k = Sgh (|3.6| ). 
With this ansatz an obvious choice for the wave function renormalization in the 
effective infrared cutoff is ZA,k = Zp^k and Zgh.k = 1- We note that the trunca- 
tion (|4.1|) corresponds to a (renormahzed) gauge fixing parameter = 1. This is 
convenient for the illustrative purpose of this section since the classical solution is 
particularly simple. For practical computations we remain more general and suggest 
the choice = 0. 

With the truncation ( [4.1|) the field equation (p.lO|) reads 

ZpFr,, + gZ'J'Z^ijj'^T^ij - ZFd^rk{-d^)A^ - Zpd^d^A: = (4.4) 

Here is a dimensionless function of {—d'^/k'^) reflecting the details of the infrared 
cutoff {V = —d'^Zp) and we choose in momentum space (cf. ( |3.12|) ) 

l-exp(-|jj 

Translating to momentum space one obtains (for details see appendix D), with 

~ 1/2 

g = gkZp^^ , the field equation 

~~^Tp I (1^^^^)^'^^^^^ + ^) = (4-6) 



d'^p d^p' 
(2^(2^ 



(2vr)^ 



+ Y(f^.''fy^u + f.^Vy.u) J (t^4 + P - P')^:i-P)A'"''iP . 



P(q) = q^ + = (4.7) 

1 -exp(-f2) 



where 



If we are interested in a polynomial expansion of Tk [i^] we can solve (|4.6| ) iteratively. 
In lowest order one finds 

{At\q)t = -~9^^ J ^,mP-\q)l'TMP + q) (4.8) 

The next term will involve a correction ~ ip'^ and so forth. (The term ~ ip^ can be 
found in appendix D (D.15)). Inserting the lowest order classical solution ( |4.8| ) into 
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the action ( [4 .11 ) and including the infrared cutoff S the effective action for the 
quarks reads 

= / {^^^(9)7^9.V^(?) + AM} + rl^)[V^] (4.9) 

with £ the Fourier transform of Ck and 

x(pi,P2,Ps,P4) = {*:(-pi)7''(r')/i;?(-P3)} {it(Plhl.(T,)M(P2)} (4.11) 

The curled brackets indicate contractions over not explicitly written indices (here 
spinor indices), i,j,k,£ = l...Nc are the colour indices and a,b = 1---Nf the flavour 
indices of the quarks. By an appropriate Fierz transformation and using the identity 

(TnAT.),' = Is'A - (4.12) 

we can split Ai into three terms [§] 

M = M^ + Mp + Mp (4.13) 
= -^{^:(-Pi)^t(P2)}{V^^(P4)^J(-P3)} 

+\ {^p:{-Pih,i'^pKp2)} {4ip^)ri'i^]i-P3)} (4.15) 

M, = -^{v5^(-pi)7,<(-p3)}{^^(p4)7'^V'?(p2)} (4.16) 
In terms of the Lorentz invariants 

s = (p^ +p2)^ = (P3 +P4)^ 

t = {Pl-P3y = {P2-P4? (4.17) 

we recognize that the quantum numbers of the fermion bilinears in A^o- correspond 
to colour singlet, flavour non-singlet scalars in the s-channel and similarly for spin- 
one mesons for A4p. In analogy to ref. we associate these terms with the scalar 
mesons of the linear cr-model and with the p-mesons. The bilinears in the last term 
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Aip correspond to a colour and flavour singlet spin-one boson in the t-channel. These 
are the quantum numbers of the pomeron. 

The effective action rfcQ[?/'] (4.9) evaluated at some conveniently chosen scale 
ko constitutes the "initial value" for the subsequent evolution of [■?/'] for k < k^. 
The evolution of Tk[ip] with the scale k is then described by the exact evolution 
equation ( p. 25 ). The latter has to be evaluated with the approximations (4J.). We 
are interested in the evolution of the two- and four-point functions for the quarks. 
The respective flow equations for these quantities obtain by taking the second and 
fourth functional derivative of eq. ( p.23| ) at ip = ip = 0. We label the different 
contributions on the r.h.s. of eq. (|3.6|) by 



d 

— rfc[^/'] = -7^-F7Av. + 7A + 7c-e (4.18) 

and discuss them separately. 
The first term 

7* = rr{(rfM + fil*T'|fir') (4.19) 

is the standard contribution of a pure fermionic theory. It is graphically represented 
by a fermion loop in fig. 1. Including only a four-fermion interaction in T^ltp] 
the contribution of 7^ to the evolution of the two- and four-point functions of the 
quarks is shown in fig. 2a and fig. 2b, respectively. In appendix E we expand the 
'(/'-dependent propagator (F^^-* [■?/'] + R^^^)~^ in the background quark fields. From 
there the contributions to the two- and four-point functions are easily computed. 
In this section we will not elaborate further on this term and rather concentrate on 
the new "correction terms" which arise from the gluon and ghost fluctuations. More 
details on 7^ can be found in appendix E (cf. (E.14) and (E.17)). 

The remaining terms '^Aipi 7a and 7c involve R^^'' and reflect the contributions 
from gluons, whereas e gives the ghost contribution. The term 



7a 



irr{(rf + (4.20) 



involves only a trace over gluonic degrees of freedom and accounts for the contribu- 
tion of gluon fluctuations around the dependent classical solution (fig. 3). The 
contribution of 7^ to the fermionic two- and four-point functions is given by the 
dependence of F^ ^^[7/^], 0]"„/ on Relevant contributions to 7^ therefore arise 
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from terms in ['?/', A, 0] which are either quadratic in A and also depend on ip or 
are cubic or higher-order in A. In our truncation the first sort of terms is absent 
and no relevant contribution to 'Ja would be present for an abelian gauge theory. 
For nonabelian gauge theories we get contributions from the three- and four-gluon 
vertices in rfc[?/',74], i.e. 

{rf%A,A^ 0))»' = Z,{{Vr{A])';^ + {vM)Z " (l>,|^l)^:} - .') (4.21) 

with 

{Vt[A])11 = -{D'[A])y,6^, + 2tg {T^)\ (4.22) 

(v,[A]y' = -{DM])D'^[A])y^ (4.23) 



Here represents the covariant derivative in the adjoint representation with 

gauge coupling g and F^'" is the nonabelian field strength associated to the gauge 
field A. We observe that A, 0] is invariant under global gauge transformations 
of ip and A. The expression for ja does not explicitly depend on ip in our truncation 
and 7a[^ = 0, A^] or rfc[0,y4fc,0] cannot contain a term linear in Ak- There is 
therefore no contribution from ja to the flow equation of the fermionic two-point 
function. An estimate of the contribution to the four-quark interaction from ja 
therefore amounts to a computation of the gluon contribution to the evolution of 
the term quadratic in A in Tk[ip = 0, A, A = 0]. This is presented in appendix F, and 
the relevant graphs are indicated in figs. 4a, 4b. Similarly the ghost contribution 
e (fig. 5) is (with the approximation ( 3.18| )) only a functional of A. The ghost 
contribution quadratic in A is also evaluated in appendix F, and its contribution 
~ ip'^ is depicted in fig. 6. 

We have computed 7^ — e in the truncation ( [4.11 ). Neglecting a contribution 
with a different fermionic index structure (the term ~ ^if^(g) in appendix F) and 
another contribution ~ ^ In Zp from the wave function renormalization Zp in Rk 
one obtains (F.56) 



GiPl -P3)M{pi,P2,P3,P4) 
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The gluon and ghost fluctuations induce therefore a running of the coefficient func- 
tion multiplying the term ~ in ( [4.9|) . More precisely, we may parametrize 

TuAi^] = -Zl J ^•••^(2^)'^(Pi+P2-P3-P4) 
{Xm{Pi,P2,P3,P4:)M{pi,p2,P3,P4) + other ludex structures} (4.25) 

Then 7a — e induces for k < ko a. deviation of from its "classical value" 

A^^nPi,P2,P3,P4) = \gl {P-\P1 -ps) + P-\p2-Pa)) (4.26) 

as given by the contribution to the flow equation 

^A„[7A - e] = ~N,gt {p-\pi - ^3)6^(^1 - Ps) + P-\p2 - Pa)<3{P2 - Pa)] 

(4.27) 

Therefore the gluon contribution to depends only on t and vanishes for t ^ k"^. 
Only the behaviour of the four-point function for t of the order of k"^ or smaller gets 
modified. Details of the behaviour of Q{q) for <^ k'^ can be found in appendix F. 

We also observe that 7^ and e account for the gluon and ghost contributions 
to the effective gluon propagator, whereas the contribution from quark loops is 
implicitly contained in 7^ (fig. 2b). We note that the latter is not distinguished any 
more from any other fermionic contributions, as, for example, from an explicit four- 
quark interaction in Ck[ip]- Only in lowest order in standard perturbation theory 
the contribution from fig. 2b corresponds exactly to the quark contribution to the 
renormalized gluon propagator. 

The contribution 7c is quadratic in the classical solution v4fc[?/'] and gives a con- 
tribution to the four-point function 

7c=^9lZl J -^^...-^^i2TT)^6{pi+P2-p3-p4)K{pi-ps)M{puP2,P3,P4) 

(4.28) 

with 

K{q) = q^P-\q) {j^r,{q) - fiFr,{q)^ (4.29) 

and 

flF = -^\riZp (4.30) 
ot 

This again results in a contribution to the running of A^ 

^A^[7c] = -\gl {K{p, - p3) + K{p2 - Pa)) (4.31) 
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Comparison with eq. (4.26) shows that this contribution accounts exactly for the 
/c-dependence of the infrared cutoff Rk contained in P. Taking only this contribution 
into acount the solution of the flow equation therefore describes the /c-dependence 
of P^^i whereas gf,^ is not modified at this level. The effective running of g1 arises 
through the contributions —7^ + 7a — e + lAip- 

The last contribution 7^^ finally involves the explicit if) dependence of the clas- 
sical solution in form of the term ~ dA^/dip. It is graphically represented in fig. 7 
and reads 

d^q ~ 2 ( ^ 



7A, = -/^^pg^U^r.(g)-r)^r,(g) 



(27r 



(27r)M2vr)4 d^J^v) Vrl'^M + i??V,, 
In lowest order we insert (f4.8|) for the classical solution 



-9 



Jip'\q)Y{T,)Mip'-q) (4.33) 



and we expand 



vf\^\ = vfn^vfm^^- (4.34) 



(rf[0])J = Z48lbl{2^fb{q-^) 

+ (4'^[0])^(g)5K2vr)^5(g-g') (4.35) 

(Here [0] is a possible contribution from C.k\^\ in eq. ( f4.1| ), as for example a quark 
mass term.) The precise index structure of r^'^''[0]'?/'?/' can be found in the appendix 
E (cf. E.4). Expanding i^f + RkY^ and assuming that L^"^ [0] is flavour- diagonal 

(4^)[0])^(g) = Z^m,(g)5,V (4.36) 

one obtains the following contribution to the evolution of the two-point function 
(fig. 7): 

(2) K-^ r, 2 / A / A r ^ . ^ ^ ~ , ^ \ 



{2tiYJ (27r)4 

'l + rf{q))-2mMl\ ,a, 
q^l + rf\q)Y + ml{q) 



P {q + p)2^Mp) , ^(^) , — ^7— ip) (4.37) 
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Two similar contributions to the four-point function arise if we include the terms 
~ ipipip in dAk/dip and similar in dAk/dip (fig. 8a). Another contribution to the 
four-point function arises through the expansion of T), [ip] (4.34) as shown in fig. 8b. 
We observe that in perturbation theory the graph 8a corresponds to a contribution 
to the renormalization of the ipipAk vertex (fig. 9a) and similarly the correspondence 
of the graph 8b is given in fig. 9b. 

To summarize, all the contributions to the flow equation (4.18) have a simple 
well-defined meaning, which can easily be expressed in terms of graphs with rules 
given in fig. 10. For small gauge coupling there is a one to one correspondence with 
associated pieces in standard perturbation theory. Our flow equation is, however, 
not limited to the perturbative regime. 

5 Heavy quark approximation 

In the limit of infinitely large quark masses our formalism simplifies considerably. If 
the momenta in the n-point functions remain bounded (and therefore much smaller 
than the heavy quark mass), we can omit in eq. (|3.23| ) the terms involving the 
inverse fermion propagator (F^^ -'[■?/'] -|- R^^^)~'^. Their contribution is suppressed by 
inverse powers of the quark masses. In the language of the last section this results 
fo 7^ = 0,7^^ = 0. The remaining computation amounts to an investigation of the 
pure gluon theory with static quarks. This is done most easily in the language where 
the gluon fields are kept explicitly and the relevant effective action is A,A = 0]. 
It is instructive, however, to understand at every step the exact equivalence with 
the effective quark theory developed in the last sections. 

If one wants to extract the effective four-quark interaction which encodes the ef- 
fective heavy quark potential, one needs the fc-dependent effective gluon propagator 
and the effective vertex ipipA. We first describe here (for arbitrary quark masses) 
the general framework how an effective four-quark interaction obtains from "gluon 
exchange" in the formulation where both quark and gluon degrees of freedom are 
kept explicitly. We then specialize to the heavy quark limit and show the equiva- 
lence with the formulation in terms of only quark degrees of freedom. The inverse 
gluon propagator 

(rf [V; = 0, A = 0, A = 0])'[ (g) = {GA{qWu + HA{q)q.qn5l (5.1) 
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specifies the term quadratic in A 

.(A) _ 1 /• A , , (2) 



= 2 / 7^^.^(-^)(r?^[^ = 0, A = 0, A = 0]):^^A^(g) (5.2) 



{2tt 

whereas the quark-gluon vertex is encoded in 



Knowledge of Ga, -f^A and permits to compute the classical solution in order 

'\\)'\\) 

{Af{q)): = -S;{q) I ^^:(p)G^(p,g)7^(T,)/V,"(p + g) (5.4) 
Here S = [T^^O] + Rk)'^ is the gluon propagator in presence of the infrared cutoff 

{RI^'^Y/Aq) = {RMK + Rk{q)quqn5l (5.5) 

and reads 

S;{q) = {GA{q) + Rk{q)r' {S^ - q'q,{HA{q) + Mq))- 

[GA{q) + Rk{q) + q\HA{q) + Rk{q))V} (5.6) 

In general, the reduced three-point function G^{p, q) may still involve Dirac matrices 
(e.g. terms ~ '^^Pu)- In the heavy quark limit considered here we can neglect this 
possibility and treat G^ as a scalar function. Inserting the classical solution into 
( ^.21) and ( ^.31) and accounting for the term A^f^^'S' ( |3.22D yields the effective quark 
four point function 



(27r)4 (27r)4 (27r 
1 r d'^pi d'^pi 



mp)r{T^\'rAp + ?)}{A'b')7.(r^).Vf (p' - g)} 



(27r)^5(pi +P2 -P3-P4) 



with 



(27r)4'"(27r)^ 

{i^l(Pl,P2,P3,P4)-A/i(Pl,P2,P3,P4) + i^2(Pl,P2,P3,P4)Ar(pi,p2,P3,P4)} (5.7) 



i^i = G^{-pi,pi-P3)G^{p4,P2-P4){Ga{pi-P3) + RkiPi-P3))~^ 
F2 = G^{-pi,pi-ps)G^{p4,p2-p4)iHAipi-P3) + RkiPi-Ps)) 
(GAiPi -Ps) + Rki.Pi - p-i)y^\GA{.Pi -Ps) + Rk{pi -Ps) 

+ (pi -p,y{HA{pi-p3)+Rk{Pl-P3)T' (5.^ 
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Here A4 is given by ( ^.131 ) and 



mp^)i^2-h)iTn,'i^tiP2)} (5.9) 

We observe that in the heavy quark approximation the coefficients of the quark 
interactions in the a, p and pomeron channel ( [4.13| ) are all given by the same function 



Fi. For the general gluon propagator discussed in this section there is also an 
additional four-quark interaction ~ A/" which vanishes in the approximation = 
Rk = used in sect. 4. The general quark bilinear is conveniently parametrized by 
the real functions Z^{q) and rha{q) 

rS = E / ^Ml)(ZMr<l>^ + ma{qmm) (5.10) 

The fc-dependence of the functions Ga-, Ha, G^, and rha relevant for the two- and 
four-point functions for the quarks can now be studied using the evolution equation 
(^) for Tk[^,A,A = 0]. 

In the truncation where only the terms ( ^.2|) , (|5.3|) and (|5.1CI|) are kept, it is 
easy to see that the contributons to the fc-dependence of G^, and rha all involve 
quark propagators. In the heavy quark limit with fixed external momenta they can 
therefore be neglected and only the fc-dependence of r^,"^-* needs to be considered. 
For Z^ and G^ we may take appropriate momentum-independent "short- distance 
couplings" 

zM = 1 

_ 1 

Gi>{p,l) = Z^{m^)g{m^) (5.11) 

with renormalized gauge coupling g taken at the scale k = and the heavy 
quark mass. We also may identify k = with the "ultraviolet cutoff" or the scale 
where the initial values for the flow equation are specified, i.e. 

Zpirn^) = 1 
GA{q;k = m^) = g2 (5.12) 

Solving the flow equation for G'a(q') for k ^ one can compute the heavy quark 
four-point function for momenta much smaller than the quark mass. For aji = 
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(see appendix F) it is fully determined by 



limG^^(pi -Pa) I^M{pi, P2, P3, Pa) + ■^^^3^^A/'(pi,p2,P3,P4)| (5.13) 

For example, a potential with a Coulomb term and a linear term with string tension 
A would correspond in ( |5.13|) to 

The flow equation for Ga is computed in appendix F. 

Let us next illustrate the equivalence with the language of sects. 3 and 4. Ne- 
glecting terms ~ {ipipY higher powers of fermion bilinears, the ^-dependence of 
Ffc ['?/'] is easily understood in the heavy quark limit: We insert in 



|r.M = |rgM + |rl5[A]|^.^, + ^F^^^fv., a],^.^, (5.15) 



k 



the identity for the classical solution 



5T, _ SA'C'S _ ^iA) 



(A). 

R)f'{q)MQ) (5.16) 



and obtain 



|r,M = ^FgM + ^rf'^'^A,] + ^F^lA,] + Ia,^-^A, (5.17) 

The first two terms on the r.h.s. of (|5.17|) do not contribute in the heavy quark limit 
with fixed external momenta, the third term corresponds exactly to 7^ — e and the 
last term to 7c. These remarks generalize, of course, to higher n-point functions for 
the quarks. 

As a side remark it should be noticed that for given vertices and propagators a 
computation of 7^1 — e + 7c as a functional of Ak is independent of the masses of 
the quarks. It amounts to a computation within the pure gluon theory without any 
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explicit reference^ to the fermions. It is therefore vahd for arbitrary quark masses. 
In consequence, the results for 7a — e reported in appendix F can be used directly 
for the r.h.s. of the evolution equation ( [4.18| ) for arbitrary quark masses. Of course, 
the expression of A^. as a functional of ip involves the fermionic part of the average 
action and therefore depends on the quark gluon couplings. Also, all vertices and 
propagators on the r.h.s. of the flow equation have to be taken for the full theory 
with light fermions. 

We conclude that the flow equation ( p. 231 ) becomes very simple if the quark 
mass is large compared to all momenta in the Green functions. This static limit is, 
however, not yet the full answer to physical questions like the decription of heavy 
quark scattering or the heavy quark potential. For these purposes the momenta 
Pi...p4 in the heavy quark four-point function have to be continued analytically to 
Minkowski space and should be taken on-shell, i.e. pi = pi = pI = pI = — 
They can therefore not be considered as small compared to m^. In particular, 
one needs in (5.8) the three-point functions G^{p,q) for p"^ = {p + g)^ = — in 
addition to Ga(5')- The flow of this vertex function can again be computed in the 
formulation with quarks and gluons. The correspondence in the pure fermionic flow 
equations is straightforward and results in a flow equation for the four-point function 
for momenta on mass shell. By Lorentz symmetry G^{p, q) can only depend on the 
three independent invariants p^, {p + qY and q^. In consequence, the on-shell vertex 
is only a function of 



Gi,ip,q)\p2=(p+q)2=_rnl = G^{q) (5.18) 
For on-shell situations we therefore should replace in ( |5.13|) g'^{rn^)G^^{q) by 
Gliq)G^^{q) = F{q) where F{q) = limfc^o^fc(g) = linifc^o ^i(g)|on sheii-The relevant 
flow equation for Fk{q) 

involves also the flow equation for G^{q). We also note that for on-shell quarks 
(^ + m^'y^)ip{p) = 0, '4'{p)(J + m^7^) = and therefore (A/") ( |5.9|) vanishes. Heavy 
quark scattering is entirely described by F{q). The heavy quark potential can be 
extracted from this function by a three-dimensional Fourier transform. 

^"^This can be generalized for a fermion field-dependent gluon propagator if terms i/jipA^ are 
included. 
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6 Conclusions 



In this paper we have developed a new formahsm of how to integrate out gluons 
in QCD. The result is an effective action for quarks with an infrared cutoff k - 
the average action for quark fields rfc[^/']. Only gluons with momenta ~ 
are included in the computation of rfc['^]. An exact non-perturbative evolution 
equation describes the dependence of F^. on the infrared cutoff k. Additional gluon 
fluctuations are included as k is lowered and lead to "residual gluon corrections" to 
the purely fermionic flow equation. Our treatment permits a smooth transition from 
a description in terms of quarks and gluons (appropriate for pcrturbative QCD) to 
effective models involving only quarks. For many practical purposes one is interested 
in the limit A; ^ where all quantum fluctuations are (formally) included. 

An approximate solution of the flow equations needs truncations of the most 
general form of F^. This restricts the present range of practical applicability of 
our formalism. With the present rough truncations in the gluon sector we see 
two areas where this method may lead to interesting new results: One is colour- 
neutral strong interaction physcis, in particular the properties of mesons. The other 
concerns heavy quark physics at intermediate distances or momentum scales, say 
((0.3 — 1) GeV)^. On the other hand, present truncations seem inappropriate 
for a quantitative treatment of confinement. Here the low momentum properties of 
the "glue" need to be understood and, most probably, additional composite degrees 
of freedom should be introduced. 

This paper is devoted to the development of the new formalism. Sect. 2 is kept 
rather general and the general field-theoretical setting can be applied to problems 
of integrating out degrees of freedom in a wide context. In the following sections 
we focus on QCD. In view of the length of this paper we have not yet addressed 
here practical applications, which will involve numerical solutions of truncated fiow 
equations. Nevertheless, the necessary analytical computations for this purpose 
have already largely been carried out and are presented in the appendices: The 
heavy quark potential involves the flow equation for the gluon propagator which 
is computed in appendix F. The residual gluon corrections for the flow equations 
for light quarks need in addition the results of appendix E. With the help of these 
formulae the way for numerical solutions of the truncated flow equations seems open 
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and we hope that interesting quantitative results for QCD at intermediate scales and 
meson physics will emerge. 



Appendix A: Field transformations in the flow equa- 
tion 

For a general fc-dependent nonlinear variable transformation = ip = 

ipk['^,ip] the flow equation (|2.9|) transforms into 

-r f(3 ^1 = -r - ^EMMM _ ^14^9^ 

dt ^^'^'^^ dt ^i*'^* 80'' dt dipi^ dt 

2 ^ ' dt ^ dt di)P dt 

dip'' dt dipl^ dt ^ ■ ^ 

(Here = <^fc[<^, V^], the partial derivative dCp'^jdt is taken for fixed and t/' and 
similar for d'lp^/dt.) The corresponding matrix of second functional derivatives f[,^'' 
reads now, with combined fields 

a^ = (y.^V'^), a^ = {r.^^). (A.2) 
.f (2)^. ^ d'^u da;, da'' dV, d-a'' 

^ ^ ' da;,da'' da; 5a' da'' da;da' ^ ' 

We can therefore write 

^r{[l,} 1^ dt )^da^' ^^-^^ 

with G the inverse of the "covariant second functional derivative" 

and "connection" 



y _ da; d'a^' da' 

"''' ~dapd^da^ ^^-^^ 



Inserting the special transformation ( p.l6| ) 



if = (pk[0, ip] = (pki-ip] + (A. 7) 
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one obtains 




(9(^° dt dif'^ dt 



The elements of the covariant second functional derivative T^^^ read 





_ d'^rk[(f,ip] 























(A.9) 



We note that eqs. (|A.8| ),( |Al9|) hold for arbitrary (fki"^] in the transformation ( |A.7| ). 
If we use in addition the definition (^ki^j] as a solution of the field equation ( p.l3| ), 
we have the properties 

Vk[^ = 0,^] = r^^] - Al^)5[^fcH] (A.IO) 

S =0 (A.ll) 

In the evolution equation for ^Vk[ip] = ^Vk[Q, + ^/s!^^^ S[ipk] the covariant second 
functional derivative reduces to the simple second functional derivative of Vk 
with respect to and ip. Also (f [.^^^ ^ and (f ^^^^^ vanish and f is therefore block- 



The ip — tp components of the second functional derivative \ V\. j are simply related 



diagonal in the and ip components, implying in turn \ V\. j ^ = (^F^ j = 



to the second functional derivative of TkUp] 
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Appendix B: Alternative formulation 



Although suggestive, the definition of [■?/'] by eqs. (|2.11|) and ( p.l3| ) is not the only 



possibility for an effective action for ip. As an alternative we propose to define ^Pkii^] 
as the solution of the field equation derived from T^l'^, ip] instead of T^lf, ip], i.e. we 
replace (|2l^) and (gl3D by 

=0 (B.l) 

Tk[^] = Tk[Mi;],i;] (B.2) 

We briefiy discuss the modifications of our formalism and compare the properties 
of Ffc ['?/'] and [■?/'] at the end of this section. The functional T^lip] is not related 
to the Legendre transform of ( |2.1CI|) , but rather to W4[J, fT] with J given 

by R^jf^ ifkl'^]. In the limit k the difference AkS between f^f^?,'?/'] and rfc[y3,?/j] 
vanishes. We conclude that Tq[i/j] and Tq[iP] coincide 

ToM = Tom (B.3) 

and both define the generating functional for the IPI Green functions for ip. Inserting 
the classical solution i^kli^] instead of (pkl'ip] in the evolution equation (A. 8), (A. 9) 
leads to three modifications: First, the last two terms in eq. (A. 8) are now absent. 
Second, the derivatives 



induce an additional contribution in (f [^^'*)^^, which now reads 



^ (2) 

Third, the off-diagonal elements of do not vanish anymore 

\^k ) - \^k ) a' g^f3 
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The role of these modifications can best be understood if we write formally the 
original evolution equation (2.9) as 

^Tk[v,i^] = ^atlndet(rl') (B.7) 

where dt acts only on k. We can express in terms of the new variables a (A. 2) 
and write 

|r4.l4a,ndet(|l(rf + «.)^^|;) (b.s) 

provided the transformation does not explicitly or implicitly involve Rk- This 
is realized for the definition of (fk by (B.l). Inserting = and using eq. (B.l), 
one obtains 

|r. w = irr {^(f f M + mr'] (B.9) 

" (2) 

in agreement with (A. 8), (B.5), (B.6). Here and Rk are matrices with a and /? 
indices, and f is block diagonal 



(f 



in'Ta = (rriv^,^])".,,.,,,,, (B.io) 

The matrix 

. (B.n) 

has diagonal and off-diagonal elements as given by 



k) ) P+ Q^. [^k ) ag^^ 



{Rk) 0-{Rk ) + ^r^[R 

[-rik) a - [-^k ) a 
{Rk) f3-{Rk ) 

{Rkf. = ^(4^^)"'. (B-12) 

Similar as to the evolution equation for Tklip] ( |2.17| ) we need information on the i/)- 
dependence of (pk and on F^ [Lpk,'4'] in order to specify the evolution equation (B.9) 
for Ffc[^/']. It is convenient to split F^ + -Rfe into a block-diagonal and off-diagonal 
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piece 



(A.)^'. = ^(R'^Y. (B.13) 

and to expand in A, using iteratively the identity 

{D + A)-^ = D-^ - D-^A{D + A)-^ (B.14) 

Even powers of A contribute only to the block- diagonal piece of {D + A)~^, i.e. 

H-^ = D-^ + D-^AD-^AD-^ + D-^AD-^AD-^AD-^AD-^ + ... 

= {D-AD-^A)-^ (B.15) 

whereas odd powers of A give an off-diagonal contribution. 

In summary, we can write the evolution equation (B.9) in the explicit form 



1 d 



iRkTAD-'f .{RkfAH-''"-' 



with 



H'^, = D%, - {RkfAD-'T'a{Rk)y (B.17) 



and 



D^ = {^fm%' + {Rkf. (B.18) 



The main difference with the formulation in sect. 2 is the absence of an effective 
infrared cutoff in the classical solution (p.lj ) if is a massless field. This solution and 
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therefore the flow equation ( [B.16| ) typically involves then nonlocalities of the form 
The effective action Tk (B.2) equals in lowest (classical) approximation the 
classical approximation to Tq and the scale k appears only in the contributions from 
quantum fluctuations. In contrast to non-local n-point functions are therefore 
already present for A; > 0. At first sight this may seem as an advantage of the 
use of Ffc since in the classical approximation the nonlocalities in have to build 
up in the course of the scale evolution instead of being included from the outset. 
Going beyond the classical approximation there is, however, a considerable price to 
pay: First, the flow equation ( [B.16| ) is more complicated than ( p.l7| ). Also at every 
scale k one has to deal with the nonlocality of the classical solution instead of the 
essentially local behaviour of F^. For these reasons we will use in the present paper 
the version ( p. 17] ) of the flow equation. 



Appendix C: Field transformation for fermionic mod- 
els 

The formalism of Appendix A can easily be generalized for fermions, with a lit- 
tle care for minus signs from the anticommutation of Grassmann variables. The 
generalization of ( [A.l| ) reads for the special case ( |A.7| ), (|A.10| ), ( |A.11| ) 



^r.w = V((ff)-^U|Ar's,„ (C.l) 



d_ 

Here the supertrace STr contains an additional minus sign for the fermionic part 
of the trace. Using STr[AB) = STr{BA) we can write (with the notations of 
appendix A) 

STr\(m-' ^] = STr [c-^^] (C.2) 



with 



and 



dt \ dt 



Here we have introduced for later convenience a factor k = with the number of 
mixed bosonic/fermionic derivatives of the type dfp* /dip or dip/ dip etc. Employing 
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(A. 11) one has 



a 

where 

is block diagonal in the bosonic and fermionic subspaces and therefore a bosonic 
quantity. The factor k, arises from the commutation of {da*/da^) with (d/da^) and 
obeys 

-1 for a* = ^*,a = ^,a* = ip,a = 
K = { (C.7) 
1 else 



where is bosonic and ip fermionic. In consequence one has k = k, and we observe 

^(2 
k 



C = f This yields finally 



where the bosonic part of Rk is simply given by R^^^ whereas the fermionic part 
reads 

= 4"^-^4"% (C.9) 

" " dip " dip 
The present formulation was adapted for Majorana spinors where ip and ip are not 
independent. For Dirac spinors the factor 1/2 in front of the fermionic trace is 
absent. This amounts to a multiplication of R^^^ by two. Furthermore, ip and ip are 
now independent and ( p.9| ) must be supplemented by a second term where the role 
of Ip and Ip is exchanged. These two identical contributions also multiply effectively 
the second term in (|C.9|) by a factor of two. The overall modification for Dirac 
spinors is a multiplication of R''^'' by a factor two. 



Appendix D: Classical field equations and propa- 
gators 

In this appendix we derive the field equations and the gauge field propagator from 
the effective average action Tk[ip,A,A = 0] in the truncation given by (^?T|). We 
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start with the gauge-invariant kinetic term 



~ 1 



F-;, = d,Al-d,Al + ~gf^/AlAl , ~g = g^Z^^ (D.2) 
and write it as a sum of terms quadratic, cubic and quartic in A 

^ kin — ^ kin,2 ^ ^ kin,3 ^ ^ kinA K'-'-'^) 



with 



r2,2 = -f J d'xA'^id.d^ - d'5i:)Al (D.4) 

^kin,4 = 9 -^C'llc J d xA'^A'^A'^A"^ (D.6) 



and 

cil = = -(T")^,(r„)^^ (D.7) 

(For the last equations we use /^^" = /"^^ = For our purposes it is 

convenient to work in momentum space where 

r(^) -^Ff z [ d% d% d% .d.f N 

r..n,3 = ^9^Ly 7 (^(^(^(27r)5(p,-p,-p3) 

^-(-pO^^(P3)^^(-P2)A^(P4). (D.IO) 
Similarly the gauge invariant fermionic kinetic term reads 
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The field equation obtains as 
6_ 



i2n) 

+f^C^.y^J -0^,0j-,^''('l + P-P')^:(-P)^'"'(P) = (D-12) 
where we use the invariant tensor 

Czxyw = fzx fywu + fzw fyxu (D.13) 

For an expansion in powers of ipil! we can solve the field equation ( p.l2| ) iteratively 
with the result 

iA'^k\qm = -9j;P''(^) J ^Mp)^'(Tz).-^r,{p + q) (D.14) 



{4\q)t = ~^9P-\q)fzy'" j 
{{q''+Pn{Af{p))l-p^{Af{p))l} 



+ {Pi-P3 - 2qrm-pihuT^iji-p3)}mP4)rT^i^{p2)}] (D.IS) 

The effective inverse gauge field propagator {r^^^)l'^{q' , q)+Rlf\q)6f6^(2'KY6{q — 
q') is related to the second functional derivatives 

52 



\A=Q 



5A^^{-q)5Al{q') 

= ZFP{q)5WA'27iY5{q-q') (D.16) 
'^^"•^ -i~gZFfrm-q)uA':,{q~q') 



6Al{-q)6Al{q') 

+ (2g' - qYA^,{q - q') - [q + q\Al{q - q')5^^} (D.17) 
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6Af^{-q)6A'^{q') 2 J {2nr 

{CJJi^A^AQ - q' + p)A-''{-p) + 2C,y,^A^^{q - + p)^-(-p)} (D.18) 

It is now straightforward to expand the effective propagator {T^^^ +R^if^)^^ in powers 
of A: 

-i~gZ^'P{q')-'P{q)-'fr- 

{{2q - q'%A>:^{q - q') + (2g' - qrA^^q - q') - {q + q\Al{q - q')Si:} 
-\fZ-/Piq'rPiqrJ^,- (D.19) 
{CJJi:Al{q - q'+p)A--{-p) + 2C,y,^A^^{q - q'+p)A:{-p)} 

-fZ-,'P{qr'P{q)-' J 0^,P(p)-'fyrf/-'. 

{{2p - - q') + i2q' - p fA^^ip - q') - {p + q\Alip - • 

{{2q - p)M^,{q -p) + {2p - qYA^,,,{q - p) - {q + p)rAl,{q - pjc^^^,} + 0(^3) 

Eqs. (D.14), (D.15), and (D.19) define the quantities needed for 7^ and 7^^ in sect. 
4. For easy reference, we also give the vertices 

'"'^ - igZFfzyw{2n)'^S{qi + q2 + qs) 



6A^'^{ql)SA''y{q2)SA'''"{qs) 

{{Qi - (l2)a5nv + (?2 - qz)^i5ua + (?3 " (lljv^i^a} (D.20) 

— , r = fZp {CJXA''''{p + q-q') 

5Af^{-q)5A-{q')5A^{-p) ^ \ u H H) 

+Cy^J^"Al{p + g - g') +Cy,JlA^^{p + q-ci)] (D.21) 



-fM27:Y5{q + p-c^-p') 



6Al{-q)6Al{q')5A^^{-p)6Al{p') 
{CJ^Si^S^, + CYJ'S^S,, + Crjl^^] (D.22) 
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Appendix E: Expansion of the effective action and 
propagator in fermion bilinears 

The effective action for quarks [■?/'] contains an equal number of ip and ip fields as a 
consequence of baryon number conservation. It can therefore be expanded in terms 
of IPI functions for an even number of fields. Denoting by a.., 7..(7..) the flavour, 
colour, and spinor indices of the Dirac spinors one has in a momentum basis 

+7 / ^■■■^,i^'2{-PMUpM\pM8{-Pz) 



A J (27r)^"'(27r)'^''" ' ^''"^^^^ 
(^^'^[OOIIS (^i'^2,P3,P4) + mi^f] (E.l) 
As a consequence of momentum conservation and global colour symmetry we write 



(r?^[0])I"(^,^') = (^)^^(2vr)^5(g-g') (E.2) 



^^^Z (^-^^'^3,P4) = (G%^^ {PuP.Ps,P.){2nr6ip,+P.-Ps-p.) 

(E.3) 



with 



-{G^%ll%iPl,-P3,-P2,PA) 

iG^)tll^P^^-P^^-P^^-P^ (E-4) 



In the approximation (|E.1| ) the second functional derivatives read 



+ 



5^T 



i^7{q + q'-p)rc{p) 
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The effective propagator (including the infrared cutoff R^^^) can now be expanded 
correspondingly. We will use here the ansatz^^ for the two-point function 

{G^'^TM = Z^{cM4+rna{q)i)5t (E.8) 

Noting that in the order {ipipY o'^^ has to take into account also the off-diagonal 
{tpip) and {tpi^) matrix elements in T\. one finds in this order: 

617 

hf''Wi:f{p)i>Up')Ci-')%<l -9+p'-p) 

'"(9 - Pi - P2) (G'''>):^^^(-<;' + P3 + P4, P3, P4) 

4 V gaimrie 

h^^'{p2 -P3-q) {G^^^J- iPi, Q+Pi- P4,P4)}/if ^(g') 

(-Pi)V'f5(P2)<(p4)V','^(-P3) (27r)'^5(g -q'-p,-p,+p,+ p^) (E.9) 



where the cutoff-propagator h reads (r^ = r^^'* in this appendix) 

" (c„(?)+r,(?))V + m^(5) ^ • 

This yields the following quadratic and quartic pieces in 7^ 



and 

(4) ^-2 f d'^PA d'^q sr^ rT(a)r,r \Uie)er \ 



^^The matrix 7 generalizes 7^ to arbitrary d. 
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(^^^XTkSi-'i^ ? - Pl - P2, -P2, -Pl){G^%iZtA-(l + P3 + P4, g,P3,P4) 



agmpe'0 ^ 



+iArfe-f3-';)(G'*'^""" 



.. .. -P3, -P2,P2 -P3 -g) 



aekpS-ff 



i^f{-pMt8{P2Wf{PAW,^{-p^){'i7:Y5{Pl +P2-PS-P4) 



where we use the shorthand (77^ = —-^InZ^) 



-2 



(E.12) 



((ca(?) + rk{q)fq^ - ml{q))i + 2m„(?)(c„(g) + rfe(?))?'7l (E.13) 



To be more concrete, we consider the truncation 



-zlj 



^{27ry S{pi + p2 - P3 - pa) 



^^J (27r)'^'"(27r)' 

{Xa{Pl,P2,P3,P4)Ma + >^p{PuP2,PZ,Pa)M p + 
>^p{Pl,P2,P3,P4)Mp + Xn{Pl,P2,P3,P4W} 



(E.14) 



with Ao.(-p4, -p3, -p2, -Pi) = Ao.(pi,p2,P3,P4) and similar for A^, Xp, X„. This yields 

Z^'(G^)i:^ip.,P2,PS,P4) = 

XM,P2,P3,PdmHSi {5'^5l - (7)^-^(7)/) 

-A.(pi, -PS, -P2,P4)mslSi (<5^<5J - (7)/(7)/) 
-^A,(pi,p2,P3,P4) ((7M)-y"(7'^)/ + (7.7)^'^(7''7)/) 

+^A,(pi, -PS, -p2,P4)S^MSi ((7.),' (7'^)/ + (7m7)/(7^7)/) 
-YXp{p„p„ps,P4)S^ASlSi{j,)^\r)s' 



+ ^Ap(pi, -P3, -P2,P4) ^^^4(7^)7^(7 



+2Xn{p,,P2,P3,P4)SA{T.),\n,'i^^ -k)s' 

+2A„(pi, -P3, -p2,P4)<^,"5S(7;)/(r^),'(?^i +1^2)i"(j^3 +1^4)/ (E.15) 
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From there we obtain 



2 ''J (27r)'^ (27r)'^ \dt 

a 

[i'l{p){-'^^^^Ncma{q){Ca{q) + rk{q))q^lK{-q, q, -p,p) 
-2^N^[{Ca{q) + rk{q)fq^ - ml{q)]<iXp{-q,q, -p,p) 



{2I - 2) [{cM + rkiq))y - mliq)) 4 (E.16) 



2 r 

+22+^m„(g)(c«(g) + rk{q))q^l\\p{-q, q - p,p) 



+2ma{q){ca{q) + rfc(g))g^(g - p)^7] A„(-g, g, 
+ [(ca(g) + rkiq)Yq^ - ml{q)]^l{p)(i^', (p) 

{-4A^(-g,p, -q,p) + 2(2 - 2^)Xp{-q,p, -q,p) + 2i+Up(-g,p, 

For an investigation of '~)f we first consider the simple truncation (4.25) which 
is relevant for sect. 4. This corresponds to 

(^^^))::::(..p.P3,..)^2zj 

{A™(pi,P2,P3,P4)5f5.^(T^)/(T,),^(7,)^^(7")/ 

-A^(pi, -P3, -P2,P4)5^5e'(T^)/m).'(7,),"(7")/} (E.17) 

We also specialize to massless quarks and Ca{q) = c{q). Furthermore, we will omit 
from now on the second term in the bracket in ( [E.12| ). (The index structure of this 
term is different and its contribution to can be computed in a similar manner.) 
One obtains 

7. = (2^- {2.Y i2nr ^ ^^^^ + P2 - P3 - P.) 

2^+^NFXmi-q, -P3, -q+Pl -P3, -pi)Xmi-q - P2 + P4,P2, -q,P4:)H{q) 

h{q - Pi + P3){2qpq^ - {pi - Pa)^?^^ - (Pi -Pa);.^^ + ((Pi?) - (Psg) - q'^)S^,^) 

{i'{-Pi)rTM-P3)Wip^)YT'i^iP2)} 

-4Am(-g,g -pi - p2, -p2, ~Pi)Xmi-q + P3 + Pi, q,P3,Pi) 

{ip{-Pih^h{q - Pi -P2){4-A -^2huTzTyiij{-p3)} 
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{i,{p,)YH{q)h,TyT^i;{p2)] 

-AX^{-q,q - pi +P3,P3, -Pi)Am(-g -P2 +P4,P2, -Q^Pa) 
{4,(-p,)Yh{q-pi+ps)i4-A +]^sHH{q)4^,TyT^Tyiji-p,)} 
mp,)rTMP2)} 

-4A^(-g, -P3, -q + Pl- P3, -pi)Xmi-q -P2+ P4, q, -P2, P4) 

{4,(-p,)^'^TM-P3)} 

mp,)YH{q)HMQ-Pi+P3W-A+AhuTyT^Tyi;ip2)}] (E.18) 

where Np is the number of quark flavours and 

H{q) = H{q)l h{q) = h{q)4 
h'\(l) = {c{q) + rk{q))q^ 

H{q) = q'h\q) - v^n{q)j = -dMo) (E.19) 

For the more general ansatz (E.12) the contribution to the flow equation of the 
four-point function is rather lengthy and may be split into different terms 

cfipi d'^pi cfiq 



{^0-0- + ^ap + ^ap + ^an + ^pp + ^pp + ^pn + ^pp + ^pn + ^nn} (E.20) 

with Acrp ~ AcrAp etc. As an example we give (again omitting the second term in 
the bracket in ( [En^ ) 

A„„ = NcXa{-q,q - Pi - P2, -P2, -pi)Xa{-q + Ps + P4,q,P3,P4.) 
E[{^a(-Pi)^!'(P2)}{^^(p4)^,n-P3)}tr(/^('^)(g-pi-p2)//('Hg)) 

a,b 

Hi^:i-pi)ii^Kp2)}mp,)ii^]{-P3)}H^^^ 

-{V'a(-Pi)^'(P2)}{V'^(p4)7^,n-P3)}tr(/i(")(g - Pi - p2)7i^^'Hg)) 
-{<(-Pi)7^-(P2)}{V^^M^,n-p3)}tr(/i(")(g-pi-j92)ff('n9)7) 

~Ki-q,P2, -q + pi +P2, -pi)Ki-q + P3 + P4, -P3, -q,P4) 

J2[{^l{-Pi)h'^''\q-Pi-P2)^':{-P3)} {i^l{P4)H'^''\q)^-{P^)} 

c 

+ {i^:i-Pi)ih^''\q-Pi -P2)ii^1i-P3)} {i^liP4)iH^''\q)i^]iP2)} 

- {^:{-Pi)h^''\q - Pi - P2)l^'i{-P3)} {4{P4)lH^''\qWP2)} 
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-Xai-q, q-pi- P2, -P2, -pi)K{,-q + P3 + P4, -P3, -g, -Pi) 

a 

+ {^l{-Pl)lr^{P2)] {^l{P4m^'\q)lh^''\q -Pi- P2)l^',{-P3)} 

- {i^:i-Pi)€iP2)} {i^l{PA)lH^''\q)h^''\q-Pi-P2W,{-P3)} 

- {^:(-Pi)7C(P2)} {4{PA)H^''\q)lh^''\q-Pi-P2)i^]{-P3)}] 

-K{,-q,P2, -q+Pl +P2, -Pl)A<x(-g +P3 +P4,g,P3,P4) 
b 

+ {V^:(-Pl)7/^^'n9 - Pi - P2)7^('ng)7V^r(P2)} {V^^'(P4)7V^)(-P3)} 

- {V'^I-PO/^^'H? - Pi - P2)7i^^'n?)^r(P2)} {^^■(P4)7V^?(-P3)} 

- {i':{-Pi)lh^'Hq-Pi-P2)H^'\q)l^nP2)} {V^^(P4)^?(-P3)} 



(E.21) 



Here the trace tr is over Dirac spinor indices. Other terms hke A^p and App can 
be obtained by inserting 7^ in appropriate places and accounting for appropriate 
factors 1/2 or 1/4. In the chiral hmit these expressions simphfy since for nia = 
0, Ca{q) = c{q) one has h{q) = h{q)4 and H{q) = H{q)4. One finds 



A 



rk{q) -mrk{q) 



q^iq - Pi - P2Y \dt 

[c{q -pi- P2) +rk{q-pi- P2)Y^[c{q) + ^^(g)]"^ 

K{,-q,q-pi -P2, -P2, -Pl)Aa(-g +P3 +P4,g,P3,P4)-M<7(Pl,P2,P3,P4) 
-2NFXa{-q, -P3, -q+Pl -P3, -Pl)Xai-q - P2 + P4,P2, -q^Pi) 

{i'a{-pi)h{q-pi+P3)i^1{p2)} {i^l{PA)H{q)r,{-P3)} 



+ - Pi +P3)7V^-(P2)} {i'i{pA)H{q)^rj{-P3)] 



(E.22) 



where we have made for the second term a replacement p2 



— consistent with 

The first term gives a contribution to the flow of whereas the second 
term, upon performing the g-integration, contributes to the evolution of Xp. We 
observe that the omitted term in Afja arising from the second term in the bracket 
in ([E. 12|) contributes only to the flow of Xp. 
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Appendix F: Evolution equation for the gauge field 
propagator 



In this appendix we compute the evolution equation for the gauge field propagator 
in a pure Yang Mills theory without fermions. We use the gauge fixing 

S,, = ^JdMd'A;f (F.l) 

for arbitrary a. We start with the gauge field- dependent inverse propagator in 
momentum space 

+tgZFf^y^{{2q~q\A>:,{q-q') 

+ {2q' - qrA^q - q') - (g + q\Al{q - q')5^^} 

+\fZp J ^,{{fJf\u + f.^V\umAl{q-q' + p)A-'^{-p) 

+2{frU.u + fJL\)A-'{q - ^^v)A-{-v)} (F.2) 

Here P^[(i) is the general kinetic operator 

P^{q) = 9{qK + Kq)q,q^ (F.3) 
with g and h functions of q^. In the classical approximation one would have 

9{q) = q\ %) = (--!) (F.4) 

a 

Here we incorporate an infrared cutoff ~ k such that for massless gluons 

lim g{q) ~ /c^ lim h{q) -. (F.5) 

For the present calculation we keep g and h as free functions. Whereas the A- 
dependent terms in (|F.2|) can be derived from a term in the action ~ ^ZpF^^Fjf''^ , 
we consider here first the most general form of the A-independent term. 

The evolution equation for the ^-dependent effective action for gluons {T^jf'' = 
Tk[A,A = 0]) 

|rl^^ = IJ^. {|[^.(g)];:[(f?V?.^(g,g)} -e = lA-e (F.6) 
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contains contributions from the gauge-field fluctuations (7^) and the ghosts (e). 
They can be identified with corresponding contributions in the flow equation (4.18) 
for the quark-effective action in the limit of a truncation where (Xk does not 
explicitly depend on ip. Defining formally a partial derivative dt acting only on P^^ 
with 

dtP!:{q) = ~Zf^^ {RkiqW: + Mq)q.q'') (F.7) 
we can write 7^ as the derivative of a one-loop expression 

7A = ^Tr4 In fi'^ (F.8) 

The evolution of the gauge- field propagator can be extracted from the term quadratic 
in A in 7^ 



(2) 1 f d'^p 

Ta = : 



^ {Va{p)A^A-p)a;{p) + Wa{p)puA:{-p)p^a;{p)} (F.9) 



(27r)^ 

Our aim is the computation of the functions Va and Wa which depend on p^. The 
ghost contributions will be added later. 
We expand 

Trdt\n{P + A) = TrdtlnP + Trdt{AP-^) - ^Trdt{AP-^AP-^) + ... (F.IO) 

where A is the A-dependent piece in f^^'' ([F.2|) . The first term amounts to an 
irrelevant constant and the contribution linear in A from the second term vanishes 
due to f/"' = 0. We obtain more explicitly 

7?^ = I J (F.ii) 

d'^q d^p pr 



dt{Ay/,iq,p)Z^\p-Xip)^Z{p,Q)ZF\P-')M} 



with 



9{p) V dw + Km 



or 



g[q) + h[q)q^ 
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Using the identity^ 



f fzyx 
J zywj 



X 

w 



(F.14) 



this yields 

,(2) 1 



Ta 



{2nY {2'kY 



d,{g-\q)[{d-l-b{q))A;{p)A^,{-p) 



q'^q.Alip)A:i-p) 



(F.15) 



9-\q)9-\q + p) {5p' + 2{pq)+2q')Al{p)A':{-p) 



+ ((4d - Q)q^q, + {2d - 3)(g>, + p'^q,) + {d - 6)p^Pu)Al,ip)A:{-p) 
-2^-^{{q' + 2{pq)fA':{p)Al{-p) 

+ [(p' - 2(m) - q^)q^qu - (g' + ^{pq)Wp. + P^qu) + gVp.]^^(p)^r(- 

6(g) 6(g + p) r 



g2 (g + p)' 



+ +(pg)Vp.] Al{p)A'^{-p)\] 



The term ~ yl^(p)A^(— p) can be projected on contributions ~ Ga using as a 
projector (^5^ — , and similarly for contributions ~ Ha with the projec- 
tor With 



/ 7^^9tVA{p,q) 



T.'i 



C~g 



(27r)' 

^ ^ dtWA{p,q) 



we obtain 



VA{p,q) = 9 ^{q) 
1 



d - 1 - 6(g) 



{27rY 

b{q) 
d- 1 



(F.16) 



:9 ^{q)9 ^{q + p) 



+ 2(pg) + 2g2 + 



_ {pq? 

g2p2 

d- 1 



p2 



+2 



2 Kg) 

d-l g2 
(pg)3 (pg)2g2- 



{(rf - 2)(g2)2 + {Ad - 6)q\pq) + qY + (4^ - 5)(pg)' 



+ 



1 b{q)b{q + p) 



d — 1 g2 (g + p)2 



{(j>^)V-(m)V} 



(F.17) 



^^The value C = Nc obtains for the gauge group SU{Nc)- We keep in this appendix the Casimir 
operator C for general nonabelian simple gauge groups. 
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and 



WA{p,q) = g ^{qy^ 



d {pqf 



1 q^ 



1 



4d-6 / 



rf- 1 



d 



+ {Ad-6)^ + d-6 



26(g) \d-2 2_^d-2q\pq) 
-q / 



q' 



d-l 



d — 1 

5d — 6 {pqY 1 (g^)^ (i {pq^q'^ 2d (pq)^ 

d-l p2 + c^-l p2 - rf-l (p2)2 ~ - 1 (p2)2 

1 h{q)h{q + p) 



d — 1 (f- (g + p)^ 



{(pg) - g p } 



(F.18) 



We next evaluate the contribution from the ghosts given by e ( p.l7| ). With the 
approximation ( p.l8|) and A = one has 

^ = Tr I [jR^t") (rQ'D,\A\ + i?f y j (F.19) 

where R^^^^ is now a function of the normal Laplacian. We generalize this formula by 
taking into account the most general effective ghost propagator. This corresponds 
in momentum space to the truncation 



^(gh){2) 



+ i?f ) \ (g', g) = ZMq)6y{2nr6{q - q') - i~gZ,^U ^A^q' - q)q 
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An expansion of the propagator up to quadratic order in A gives 



(gh){2) 



(gh) 



{q',q) = Z:^P-^\q)5l{27^)%q-q') 



M~gZ^,lP,,}{q')P,^\q)U ^q"A-{q' - q) 

-fz^^P,^\q')P£{q)UJ.\- 

I (l^^gh'(p)^>^^M(9'-p)^^(p-?) 



This yields, up to an irrelevant constant 



Cg 



7.2 



{2t:Y dt 



Rf\q)P^'iq) 



I J^^gh'(^-P)^'(^-P)'^^^(P)^^(-P) (F.22) 



(2vr) 



(F.20) 



(F.21) 
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With 



(27r)° 



and defining g^^i and /igh in analogy to ( [1:''.16| ) one finds 



C~g 



d- 1 



p2 



d'^q ~ 



(2vr) 



~dt{P^,l{q)P^^\q + p)q''{q + pr} 



^gh = ^^^gh'(?)^gh'(9 + P) 



p2 



(F.23) 

(F.24) 
(F.25) 



gh 



C~g 



7.2 



1 1 

d — Ip"^ 



1 



(27r 



{^gh'(9)^gh'(? + + 

(F.26) 

{pqf 

--HI/"/ — II — ^ 



,2^2 



pZ pZ 



(F.27) 



«^gh = ^^^gh^(?)^gh^(? + P) 

Combining the contributions from 7^ and e we arrive at the fiow equation for 
the gluon propagator for the pure Yang-Mills theory. It describes the /c-dependence 
of the function G'a(q') and HA^q) as defined in (^.l]) and reads 

^Ga{p) = q-^{Zf9{p) - Rk{p)) = Va{p) + V,^{p) 

-iJ^b) = gliZphip) - Rk{p)) = Wa{p) + WM (F.28) 

We observe that Va and Wa are defined in terms of Ga and Ha through the relation 
T^^\0) = t^^\0) - R\^\ Ga = Zpg - Rk, Ha = Zph - Rt- The system ( ^1281) con- 
stitutes therefore a coupled system of partial differential equations for the functions 
Ga and Ha which depend on two variables, namely k and p^. It describes the scale- 
dependence of the gluon propagator under the infiuence of the three- and four-gluon 
interactions contained in {F^^f and the corresponding ghost contributions. 

It is instructive to study first a simplified version of this system where the prop- 
agator is approximated on the r.h.s. by 



Zpgiq) = ZFq^ + Rk{q) = ZpP{q), hiq) = 



(F.29) 
(F.30) 



This corresponds to the classical approximation ( |I''.4| ) wit a = 1. It implies 6(g) = 
and 



VA{p,q) 



-\p^\q)P-\q + p) 
+ {d-l)p-\q) 



6d-8 2 ^, , n 4d-6{pqf 



(F.31) 
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{M - 6) 



(pq) 



+ d-6 + 



2d(2d - 3) {pqf 4d - 6 g2 



p^ d — 1 (p^)^ d — 1 p'^ 

Combining this with the ghost contribution one obtains 

(9 _ , , , , ^ , o r d'^q d 



dt 



(F.32) 



(F.33) 



Qd - 10 
d-l 



a — 1 p"' 



d_ 
di 



d-6 + A{d-2)^-A 



pq d — 2 q"^ 



+ 4d 



- 2 {pqf 



p^ 



d-ip^ d-i{p^y 



(F.34) 



If we want to study the behaviour of the gluon propagator for small values of the 
momentum p"^ ^ k"^ we can expand the functions Ga and Ha in powers of p^. For 
this purpose we expand P~^{q + p) (with P(q) = -^P[q)) 

P-\q+p) = P-\q) - 2{pq)P{q)p-'{q) 

-p'P{q)p-\q) - 2{pqf{P{q)p-\q) - 2P\q)p-\q)) + Q{p^) (F.35) 



and use the identities 



(27r) 

d'^q 



(F.36) 



/d q 1 f d q 

j^f{q^)qi.quqpqa = ^(^^(V'^p<^ + '^mp^^-^ + V'^^^p) J j^i'^f^i^^ 



(27r)' 
This yields 



d_ 
di 



G'^(p) = T/(0) + y(V + O(p') 



(27r)^ 



with 



(27r)<=' dt 



d 



(F.37) 
(F.38) 



2^ J {2T:Ydt 



d'^q d [^„_2/ V 6(i-4 On/ x n-s/ x 
' 5P ^{q) —q^P{q)P \q) 



d 



^^^"^ %'[P{q)P-\q) - 2P\q)p-\q)] 



d{d + 2) 



(F.39) 
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and 

'''"~\[P{q)p-%q)-2P'iq)P'\q)]] (F.40) 



ci(ci + 2) 

With X = we define the constants /,^, mj^, Vd by 



1 r'^ ri r) 

id _ , 2n-<i / ^ 



= --fc^"-'^' / dx x^-'—P-"" (F.41) 
2 Jo at 



and find 



dxx2^{p2p-"} (F.42) 
V-,' = 2'^+Vir(^) (F.43) 

= -AvdCfk''-^ (^{d - - ^^^f (F.44) 



d(d + 2) 
= 2t;rfC^2A;^-^ ( {d - 6)/^ - 



-16v, 



•d- 



d{d + 2 



Cf \ dx^^ |xf+ipp-n (F.46) 
J dxot ^ ' 



,2n-2 



'd{d + 2) 

Using the definition ( |4.7| ) for P{q) one easily verifies the identities = 1, 
1. For d = 2 the coefficients V"'-'^^ and 1^^°) vanish and V^^^ is proportional to 
5ll - Ami = 10 In 2 - 2. For d = 4 one has = l/327r2 and 

V^'^ = -^^fe{3lt-2lt)=0 



167r2 



= (F.47) 

A similar expansion for small values of can be made without the approxima- 
tions ( F.29|) . One replaces P in ( p.41| ), ([F.42|) by Pa or Pgh and adds a contribution 
~ din Zp/dt for the t-derivative of the factor Zp in P^. 
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For a study of the behaviour of ( p:''.33| ), ( H:''.34D for general p it is convenient to 
introduce in the definition of P an exphcit ultraviolet cutoff A^'^ 



exp [-^) - exp (-f^) 



/ daexp{-aq^) (FAS) 



This permits to perform explicitly the momentum integration 



^GM = {d- l)Cf^^ d^j^^ exp(-ag^) 



d 



s / 

2 Ot J A 



da 



dp 



(27 

exp(— ag^ — /3(g + pY 



6d- 10 



q^ + 2{pq) + - 4 



d'^q 
{2'kY 
d-2{pq){pq) 



d-1 ^ d-l 
-2{d-l){AiT)-iC~gV~^ 



-i(4.) 



^2 



3d -A 2 

a + p 



d_ 
di 

U 



A-2 



dad(3{a + [3)' 
2a(3 



exp 



al3 
a + (3 , 



(F.49) 



(^-6-4(c^-2)(^)exp(--^/) (F.50) 



Using the integral 

a + p 



(F.51) 



with 



J(y + l) = eA^Jiy) (F.52) 



exp- 



and 



-"^•^The previous form of P{q) obtains for A ^ 00. 
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+ ::^ h^P 7:^-7 exPTTT? (F.53) 



we find 

d 



g{p) = -2(d- l)(47r)-2fc' 



.(i-2 



+2(47r)-i^'^-^ exp (^-|^) |(3ci - 4) - 1) 



-4(d - 2) J - 1) + 4(rf - 2) J I (F.55) 
In particular, this yields for d = 4 

§^GM = Pa(p) (F.56) 
- 1^*^^^ - - - (-^) + - 8^) } 

The qualitative behaviour of Pa{p) can easily be understood by expanding for <^ 

. 2/10 15 , 79 (p^y 61 /pn^ 29 /pn^ 

^^^^^ 167r2^ ^ 1^ 3 16 A;2 ^ 480 [k^ J 2880 l^A;^ j 13440 l^A;^ J ' 

(F.57) 

whereas for ^ 2k'^ one has 

It is instructive to write the gluon propagator in terms of momentum-dependent 
wave function renormalizations Z^, Zi for the transversal and longitudinal gluons 
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and a gluon mass term m\: 

Ga{p)S^,u + HA{p)p^,Pu = rh\5f,y + Zt{p){p^5f,y - p^p^) + Zi{p)pf,p^ (F.59) 

m\ = Ga{0) 

Ga{p) - G'a(O) 



Zt{p) 



p2 



ZM = HAp) + (F.60) 

We may now define a fc-dependent gauge-fixing parameter a{k) and Zp by 

1 



a{k) 

Zpik) = Zi(0) (F.61) 

where both definitions make only sense for ^z(O) and ^t(O) positive. From ([F.47|) 
we obtain 



dt a 

^Zp = = -C-^Zp 
dt 3 167r2 

= V^(°) = (F.62) 

where we should mention that the vanishing of V^^^ is particular to our choice of 
the cutoff function Rk. Even though a{k) does not depend on k in our approxima- 
tion, the relevant quantity for the validity of the approximation (|F.29|) is the ratio 
Zi{0)/Zt{0). The renormalized gauge fixing parameter 

aR{k) = ZF{k)a{k) (F.63) 

obeys 

—an = —C^^aR (F.64) 

dt 3 167r2 " ^ ' 



and moves away from = 1 appropriate for ( |t''.29| ). The approximation ( [t''.29| ) 
will therefore break down for k in the vicinity of the confinement scale where g1 
becomes large. We finally give the running of g"^ in terms of the j3 function for the 
renormalized gauge coupling g\ 
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For comparison, we have also computed the flow of the gluon propagator for 
a = 0. We start by summarizing the preceding results for arbitrary a. 



and 



where 



d_ 
dt 



GA{q) = KglZp 



(2vr) 



3 - ^b{q') 



^ (GA(g') + Rk{q')y' (GAiq + q) + Rk{q + g ))" Z 



5q +2{qq) + —q - y 

(qq'r 



/2 

q - 



q 

{qq'f 
q2 



/2 



(2vr) 



(G^(g') + Rk{q')y' {GA{q + q') + Rk{q + Z 

40{qqT 10 q'' jqq') 
3 (g2)2 3 g2 + ^2 



2 4(gg') 14 (gg')^ 1 g'^ 4 (gg') 



/\3 



+ 



-h{q')b{q + q') ^ , — , 
3 (g + g')^ \ g' 



3 g2g'2 ' 3 g2 



j jqq') 

3(g2)2 3(g2)2g'2^ 



'^P,^\q')P,^\q + q') 



h{q) 



3M+4 

g^ 



-'I 

(49? 



,/2n 



(g2)2 



(F.66) 



(F.67) 



(F.68) 



GA{q) + Rk{q) + {HA{q) + Rk{qW 
For a — > we observe that Ha diverges ~ ^ and 6(g) approaches one. Defining 
as above a renormalized one sees that = is now an infrared stable fixed 
point[^ We further observe that for ur — > the function F2 defined in eq. ( [5.8| ) 
equals Fi/{jpi — psY and only the function GA^q) determines the effective heavy 



We choose Rkq^ = ( ^7 l) Rk such that 6(g) = 1 + 0(q!_r) for all values of g^ 
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four-quark interaction. The flow equation for Ga can be written in a more compact 
form using 

PA{q) = Zp\GA{q) + Rk{q)) (F.69) 
and we approximate the ghost part of by 



(F.70) 



This yields the flow euqation for the propagator for = 



§-^GA{q) = N^glZp 



~-^PA\q')PA\q + q') 



(2vr)4- 
13g2 - 14(gg') + lOg'^ 



-10 



q" 



^ {qq'f ^ q^ q^q'"^ - 



J\2 



q 



l2 



(F.71) 



{q + q'Y 

The r.h.s. of the flow equation (|F.71D involves not only GA{q) but also the renor- 
malized gauge coupling and the gluon wave function renormalization constant Zp. 
We approximate Zp by the coefficient of the q^ term in Ga- More precisely, we ex- 
pand for small g^, GA{q) = fh\ + G^^q^ + G^^(g^)^ + and identify Zp with G^^ . 
For the flow equation for the renormalized gauge coupling we rely on the fact 
that the first two coefficients of the j3 function 



dt 



Cl 



22Nr 



-Cl 



at _ 

167r2 
204 



C2- 



9l 



IGvr 



2^2 



C2 



9 " 



(F.72) 



are universal In the region of large g^ we may also use nonperturbative estimates of 
(3g2 derived by different methods . An ansatz for f3g2 combined with an estimate of 
fjp fixes also the evolution of g"^ and provides all information needed for a numerical 
solution of the flow equation ( |1''.71| ). 



References 

[1] Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122 (1961) 345 
[2] W. Buchmiiller, Quarkonia, North Holland, 1992 



54 



[3] U. EUwanger and C. Wetterich, Nucl. Phys. B423 (1994) 137 

[4] C. Wetterich, Nucl. Phys. B352 (1991) 529, Z. Phys. C57 (1993) 451, C60 
(1993) 461 

[5] C. Wetterich, Phys. Lett B301 (1993) 90 

[6] D. Jungnickel and C. Wetterich, Heidelberg preprint HD-THEP-95-7 

[7] M. Renter and C. Wetterich, Nucl. Phys. B417 (1994) 181; B427 (1994) 291 

[8] C. Wetterich, Z. Phys. C48 (1990) 693 

[9] C. Becchi, on the construction of renormalized quantum field theories using 
renormalization group techniques, in Elementary Particles, Field Theory and 
Statistical Mechanics, eds. M. Bonini, G. Marchesini, and E. Onofri, Parma 
University, 1993; 

M. Bonini, M. D'Attanasio and G. Marchesini, Nucl. Phys. B418 (1994) 81; 
B421 (1994) 429 

[10] U. EUwanger, Phys. Lett. B335 (1994) 364; 

U. EUwanger, M. Hirsch, and A. Weber, preprint LPTHE Orsay 95-39 



55 



